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Abstract 

We emphasize that a muon based neutrino factory could show the existence 
of light right-handed neutrinos, if a deficit in the number of detected events is 
observed at a near detector. This could be as large as ~ 10% if the size of the new 
interactions saturates the present limits from electroweak precision data, what 
is not excluded by the oscillation experiments performed up to now. A simple 
model realizing such a scenario can be obtained adding right-handed neutrinos 
to the minimal Standard Model, together with an extra scalar doublet and a 
triplet of hypercharge 1. In this case, however, the possible deficit is reduced by 
a factor of ~ 3, and the Yukawa couplings must be adequately chosen. This is 
also generically required if lepton flavour violation must be below present bounds. 

1 Introduction 

Neutrino oscillations provide the only observational evidence of new physics beyond the 
minimal Standard Model (SM). At present these oscillations can be fully explained by 
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introducing neutrino masses and the corresponding charged current mixing [T]. These 
observations, however, cannot distinguish between Dirac or Majorana neutrinos, nor 
they require new interactions (NI) 0. The obvious question is then, where do we have 
to look in order to determine the neutrino character and/or to observe possible NI 
involving light neutrinos? This will become especially relevant when we face the need 
to interpret new data with higher statistics and precision, as foreseen from a neutrino 
factory [3] . In the following we show that present experimental constraints leave room, 
corresponding to a ~ 10% deficit in the expected number of events in appropriate 
processes, for observing NI involving light right-handed (RH) neutrinos. Such a sce- 
nario can be easily realised with a mild extension of the SM, through the addition of 
additional scalar weak isodoublet (to be denoted by rf) , and a scalar iso-triplet of hyper- 
charge 1, (denoted by A), besides three RH neutrinos, v % R . But it requires an adequate 
choice of Yukawa couplings to suppress lepton flavour (LF) violation; moreover, for this 
particular model, the deficit allowed by current electroweak precision data (EWPD) is 
reduced by a factor of ~ 3 compared to the general case above where arbitrary NI are 
parameterized by gauge-invariant dimension six operators with unrelated coefficients. 

Within the SM muons only decay into left-handed (LH) neutrinos. Even if the 
spectrum is enlarged to include their RH counterparts, these are not produced in such 
decays because they have no gauge interactions, and neutrino masses are negligible. 
On the other hand, if other interactions are present in nature, a muon based neutrino 
factory could inject an admixture of neutrinos with both chiralities. We will show 
below that the limits on NI involving RH neutrinos are to a large extent those derived 
from (inverse) muon decay, and therefore relatively weak. This then is a promising 
reaction where to look for new physics effects in the neutrino system. 

Let us first, however, state our setup. We assume that the three light neutrinos 
are Dirac-type neutrinos, i.e. that there are three light neutrino singlets beyond the 
minimal SM, and that lepton number (LN) is conserved. In practice this is not a re- 
striction on the light neutrino character, but on the type of NI. Notice that neutrino 
masses are negligible in all experiments performed up to now, except in neutrino oscil- 
lations (and eventually in neutrino-less double j3 decay, Oz//5/3). Thus, we can assume 
that all interactions conserve LN because neutrino masses are much smaller than the 
energy relevant in the processes considered and/or the experimental precision is much 
lower than the size of the effects proportional to them, as it is the case for all foreseen 
experiments not involving neutrino oscillations and excluding 0^/3/3. Though we could 
consider LN violating NI, their effects can be ignored in our analysis d, as can be the 

1 In what follows we will ignore the LSND data [2]. 

2 An effective theory only involving the light SM fields and invariant under the SM gauge group has 
only one dimension five operator violating LN [4], the famous Weinberg operator [5] giving Majorana 
masses to the LH neutrinos and then negligibly small. There is no operator violating baryon minus 
lepton (B— L) number of dimension six [BJ. Thus, operators of dimension six violating LN also violate 
baryon number (BN), and then involve quarks and will play no role in our analysis. If the effective 
theory also includes RH neutrinos, as in our case, there are two additional dimension 5 LN violating 
operators [7], one of them generates a magnetic coupling for the RH neutrinos and is very strongly 
constrained when the vr are light [5], the other generates a correction to the vr Majorana mass term 
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LN violating effects from neutrino masses. 

Hence, the NI effects we are interested in will be relatively large and LN conserv- 
ing; whereas neutrino masses will be safely taken to vanish. The effective Lagrangian 
describing such a scenario will not distinguish between (i.e. approximates equally well) 
the case of exact LN conservation with very small Dirac masses for the three light neu- 
trinos, and the case of negligible Majorana masses for the six light neutral fermions, 
the only vestige of the very slightly broken LN in this case. This was explicitly proved 
in [S] for (inverse) muon decay assuming no additional constraints on NI. Note that by 
similar arguments we can also neglect LF violation induced by light neutrino mixing 
(to a very good approximation). 

The most general four-fermion effective Hamiltonian describing muon decay reads 

n^ u ve = Yl sip (e r 7 v e a ) (v$ r 7 ^ + h.c. , ^ 

a, j3 — L, R 
7 = S, V, T 

where «, (3 label the chirality of the neutrinos, while 7 refers to the Lorentz character 
of the interaction (scalar, vector and tensor). The present limits on the size of the 
various coefficients will be discussed in the following section, here we merely remark 
that the two couplings gY L , g% R are also associated with the largest departure from 
the SM predictions for the number of events to be detected by a neutrino factory. 

There are many other available electroweak precision data that can be used to indi- 
rectly constrain the Hamiltonian ([1]), and in particular those two couplings. To derive 
such restrictions one can take two routes. The first one consists in re-writing ([1]) as a 
linear combination of higher-dimensional operators invariant under SU (3) c x SU (2) L x 
U (l) y , and adding these operators to the SM Lagrangian [TU]; the coefficients of the 
resulting effective theory can then be bound using experimental data at all available 
energies. In this approach the physics responsible for generating these operators is 
left unspecified, except for the requirement that its characteristic scale lie beyond the 
electroweak scale. The advantage of this approach is its generality, since it is not tied 
to any specific assumption about the physics beyond the SM, its disadvantage is the 
proliferation of coefficients, and the fact that more than one gauge-invariant operator 
contributes to each term in ([T]). 

The second route is to extend the SM by adding a specific set of new fields (such 
as the rj and A mentioend previously) and interactions. This has the advantage of 
providing a specific scenario for the physics beyond the SM, but the results obtained 
are often specific to the assumptions made in constructing the model. At scales below 
those of the heavy particles this model will reduce to an effective theory of the type 
mentioned above, except that the effective-operator coefficients are all expressed in 
terms of a small number of parameters and can be constrained more tightly. 

and is therefore also negligible. There is only one dimension six operator violating B— L (and LN) but 
involves four RH neutrinos [7], and then is uninteresting for us. The other dimension six operators 
violating LN also violate BN, and can be also ignored. 
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In the following we will examine both of these possibilities. In the next section we 
consider effective Lagrangian approach, which we denote by El, where we choose a set 
of 4-fermion effective operators that generate ([T]) at low energies and have little impact 
on other electroweak observables. In section [3] we will also consider a specific extension 
of the SM, which we denote by E2, based on an extended scalar sector. Finally, in 
the last two sections we discuss the implications of these SM extensions for neutrino 
oscillations and other experiments, respectively. 



2 Electroweak precision data constraints 

One effective-Lagrangian extension of the SM, which denote by El, consists in adding 
to the SM the following set of effective operators 

AG F 



gLm e L)(i>R)+g S RR (M) « 



y/2 

+- 9lr (^£) ™* te) - g s RL (M) iV2 fefl 



+ 5j f {1%-flQ («) + 9rr (WT%) (<l^n) (2) 



v 



9lr 

V 2 



h.c. 



where stands for one of the three SM lepton doublets (/ = e, /i, r) and <p is the 
SM Higgs doublet, and v ~ 246 GeV its vacuum expectation value. We do not use the 
basis proposed in [1] for writing the beyond the SM dimension six operators; still it is 
important to note that this basis must be extended to include the light RH neutrinos 
[7]1. El is manifestly invariant under SU (3) c x SU (2) L x U (1) Y , and LF (and LN) 
conserving. Note that g\ R and g\ L are associated with operators of dimension eight: 
VrI^^R an d ^rI^^r have hypercharges Y — — 1 and 1, respectively, and we cannot 
construct a gauge invariant vector with opposite hypercharge using only II and II] at 
low energies ([2} reduces to ([I]) with the definition g}^ = 1 + 5g\ L . 

In (j2J) we have chosen to write all dimensional coefficients in terms of v. This implies 
that the natural size for the coefficients is 

J (t>/A) 4 for > ^ 

y ^ \ (*;/A) 2 otherwise , 1 1 

where A denotes the heavy scale of the physics responsible for generating the corre- 
sponding operator. 



3 Although in that reference the RH neutrinos are assumed to have masses of few hundreds of GeV. 
(Note that the operator OqncIQ in Eq. (7) of that reference is redundant.) 
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A characteristic feature of this particular extension is that, except for the effect on 
the muon decay constant and inverse muon decay, EWPD are blind to the operators 
in Eq. (T2J). For instance, although gf R and g\ R contribute to e + e~ — > V^y^ and affect 
the Z° invisible width, the effects is negligible compared to the Z° pole contribution. 
El does not contribute to LF violating processes either, because it does not include 
LF violating operators. Similarly, universality is preserved since the gauge couplings 
to leptons remain the same as in the SM. 

It must be emphasized, however, that in specific models (such as E2, described 
in section [3] below) the various coefficients in Eq. ([1]) are written in terms of the 
fundamental parameters of the theory, these parameters appear in all other interactions, 
and in general will contribute to other observables, which in turn can be used to 
stringently constrain g^*. In fact, it is non-trivial to find models where these limits are 
not so strict they exclude further observable effects from ([I]). 

As we argue below the largest departure from the SM predictions is parame- 
terised by g\ L and g RR . Constraints on these two operators (ifyjQ (^rIl) and 

(^l7^l) yLl^l) i R Eq. © can be obtained using current data, as described in the 
Appendix. At 90% C.L. we find 



case 


\gL = l + 6gY L \ 


\9 S rr\ 


(a) 


> 0.960 


< 0.550 


(b) 


> 0.957 


< 0.579 


(c) 


> 0.9998 


< 0.054 



(4) 



where the limits in case (a) are obtained directly from using muon decay and 
z/e — > ujj, data [U [TT] ; in case (b) from a global fit using precision data to the operator 
coefficients with the SM parameters fixed at their best-fit values 0; in case (c) as in 
case (b) but taking one effective-operator coefficient to be non-zero at a time. This last 
possibility, though often adopted for simplicity is seldom realistic: as noted previously 
one must expect the heavy physics to generate several operators with related coeffi- 
cients, so that fits allowing for several non-zero interactions become compulsory. Even 
more, if the new physics is to have sizable indirect effects, then it must also conspire 
to preserve the excellent fit of the SM to the EWPD at the 0.1% level 0. In Fig. Q] we 
plot the bounds for case (b). The results in (j2J) indicate that the two parameters in the 
fit are highly correlated. This is expected since the fit is dominated by the constraint 
on the strength of the muon decay constant G^, which is proportional to the SM Fermi 

4 See |12j and reference therein for a fit to operators involving only left-handed neutrinos. 
5 Thcsc limits stay mainly unchanged if the SM parameters are also left free in the global fit. 
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Figure 1: 90% C.L. bounds for El case (b) in Eq. (j4]); and the same for the SM 
extension E2 in next section. The narrow bands between the origin and the crosses 
define the 90% confidence region for the global fit to the two parameters for El (left 
cross) and E2 (right cross) respectively. 



constant Gp in Eq. ([T]): 



G^ — Gf 
1 



r V I 2 I V I 2 I V I 2 I V I 2 

$9ll\ + \9rr\ + \9lr\ + \9rl\ 



+ ^bLr+b^r+bLr+bBLr) & 



+ l(\9 T RL 



I 2 




+ 


9rr\ 


2 


T 1 


+ 


9lr\ 



= G F A 



with the proportionality constant A restricted by the global fit to the interval 

0.9997 < A < 1.0004 at 90% C.L. . (6) 
This bound alone constrains the coefficients to a narrow band in the bg\ L — \9%j{\ plane 

\9 S rr\ 2 * -8 SgL , (7) 

as depicted in Fig. [TJ 

Using the bounds above we can derive limits on the scale of new physics responsible 
for the two operators being considered. Using ([3]) and assuming the underlying physics 
is weakly coupled we find the weak constraints A > 130 GeV for g^R and A > 500 GeV 

for 89ll- 

As stressed in the introduction, we are interested in those interactions in the muon 
decay effective Hamiltonian that could allow for the largest deviation from the SM 
prediction of the number of events that may be eventually detected at a neutrino 
factory. In general, while it is clear that a negative &g\ L is strongly favored for this to 



6 



be possible (see Eqs. (JSJ Ej)), one may wonder whether any of the other interactions 
involving a RH neutrino may play the role of g RR . This could be the case for g RR or the 
LR and RL operators in Eq. ([2]). However, unlike for g RR , which cannot be separated 
from bg\ L in muon decay experiments since we do not measure the polarization of the 
final neutrinos, all the other interactions are contrained by the absence of any significant 
deviation from the V— A prediction in the spectrum of the outgoing electrons. In 
particular, \g\ R \ < 0.034. On the other hand, the bounds for some of the LR and 
RL interactions are relatively weak, but they also generate radiative corrections to 
the neutrino masses and, assuming naturality, are constrained by the associated upper 
limits P3]: \g s LR \ , \g^ L \ , \gl R \ < 10~ 2 and |#f L | , \g% R \ , \g^ L \ < 10~ 4 . Therefore, when 
we consider the data from inverse muon decay, where the incident neutrinos come from 
pion decays and then are LH, and derive a relatively weak bound on 5gY L , this can be 
only compensated by g RR in order to satisfy ([H]), but in contrast with the remaining 
operators in Eq. ([2]) has no further constraints. 



3 A simple SM extension 



Let us discuss a simple model realising the former scenario. It extends the SM including 
besides three RH neutrinos v % R with zero LN, a second scalar doublet rj with LN equal 
to —1, and a scalar triplet A with hypercharge 1 and LN equal to —2, 

V + \ A _ ( A+ V2A+^ 
i) /' \ V2A° -A+ 



V 



(8) 



as in [13] and [32], respectively. They both can acquire a vacuum expectation value 
through very small LN violating couplings to the SM Higgs (which carries zero LN), 
and provide masses and mixings to the light neutrinos, which must be then fit the 
observed spectrum in oscillate experiments B Bnt, these couplings , as the neutrino 
masses, can be safely neglected in our analysis: we can assume that light neutrinos are 
massless, and that LN and LF are both conserved. Thus, besides the kinetic term for 
the RH neutrinos, the SM fermionic Lagrangian acquires only two more terms 



fijU R l jT ia 2 rj 



\ ij l i LCia 2 Al j L + h.c. 



(9) 



with only f ee ^ and A e/X = A Me potentially large enough to produce measurable effects 
in current and near-future experiments. The integration of the extra scalars out results, 
in particular, in the following contributions to the muon decay effective Hamiltonian: 



4G F 

with coefficients 



9rrC 1 



v R l L 



+ 



+ h.c. 



9rr 



feefw 



2 v / 2Gi 



M 2 



and 



SgL 



A 1 " A 



V2d 



Mi 



(10) 



6 There can be also tiny RH neutrino mass terms. 
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respectively, where M Vj ^ stand for the scalar masses. The full set of dimension six 
operators arising from the integration of the scalar triplet A is given in [TS] . 

In this case, named E2 in the former section, the EWPD analysis presented in the 
Appendix implies more stringent bounds on Sg^L an d gfm than for El: 

I 1 + 5g v LL \ > 0.988, and |^| < 0.313 . (12) 

The corresponding band is plotted in Fig. [TJ The bounds obtained are tighter because, 
as emphasized previously, the integration of definite new physics also gives, in general, 
operators contributing to other processes, further restricting the model. In the present 

case the integration of the A generates also the operator (/£7 M /|) (IlI^l) > w hich has 

the same coefficient &g\ L as (^f,7 M ^) ^xT/Jfj , an d contributes to v^e — > z/^e, further 

restricting the allowed deviation from the SM predictions 0. 

LF violation is below experimental bounds because similarly to the El case the only 
non-negligible couplings in Eq. are f ee ^ and A eM = A^e, and because the scalar 
doublets and triplet mix very little, as required by approximate LN conservation. The 
absence of new r couplings and that the SM gauge couplings stay unchanged guarantee 
the agreement with universality constraints on the lepton sector. 

4 Neutrino factory predictions 

The relevant phenomenological question is where could the RH neutrinos be eventu- 
ally observed if the 5gJ jL and g RR four-fermion interactions are non-zero. Obviously, 
they can be probed in a more precise inverse muon decay experiment: a more precise 
measurement of this process could give evidence for those NI (or reduce the allowed 
deviation from the SM in Eqs. fj4|TT2T) and Fig. [1]). But a muon-based neutrino factory 
will also be sensitive to them. Indeed, if a substantial amount of the neutrinos pro- 
duced in muon decay are RH, a near-detector sensitive to neutrino-hadron collisions 
will observe a deficit in the same proportion, and this deficit would be twice as large if 
the detector could also measure the inverse muon decay process. In Tabled] we give the 
maximum deficits expected in the case that the new interactions saturate the 90% CL. 
bounds obtained in the previous sections. Whereas in Fig. [2] we show the predicted 
deficit as a function of bg\ L for the SM completions considered. 

If the precision in the measurement of the inverse muon decay is improved by a factor 
of 2 without modifying the central value, the limits on the NI would be strengthened and 
the allowed deficit of observed events at a neutrino factory reduced; the corresponding 
percentages in such a case are given in parentheses in Table [TJ For the El extension the 
deficit would be reduced to ~ 5%, being the 90% CL. bounds in this case |l + 5gY L \ > 
0.975 and |^| < 0.442. In the E2 case the improvement in the precision of inverse 

7 In our approximation (implying negligible LN violation in the scalar sector) there are no tree-level 
contributions to the oblique parameters. 
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Process 


El 


E2 


v-N 


8.5% (5.0%) 


2.5% (2.5%) 


IMD 


15.4% (9.3%) 


4.8% (4.8%) 



Table 1: Maximum deficit in the number of observed events in a near detector sensitive 
to neutrino-nucleon collisions {v-N) and inverse muon decay (IMD) for the two SM 
extensions discussed in the text. In parentheses we show the deficits expected in the 
case that the precision on the measurement of inverse muon decay is improved by a 
factor of 2. 

1.00 
0.95 
| 0.90 
0.85 
0.80 

-0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 



Figure 2: Percentage of detected events in a near detector at a neutrino factory com- 
pared to the predicted number by the SM, as a function of the NI strength 5g\ L . 
The solid (dotted-dashed) curve corresponds to neutrino-nucleon (inverse muon decay) 
collisions. The vertical lines stand for the El and E2 limits on 8g\ L m t ne text. 



muon decay would have no appreciable effect because this constraint would be weaker 
than the one derived from ev elastic scattering. 

For a neutrino factory to be sensitive to a deficit < 3% (the smallest value listed 
in Tabled]), the neutrino flux must be known with enough precision. Besides, large 
fluxes are also required to have a large number of events at a near detector in order 
to keep the statistical error small. For a detailed study about future neutrino factories 
see [3]. Assuming 10 21 muon decays in one year of operation, the number of expected 
v-N events at a near detector (such as the one described in Table 1 of [3]) is of the 
order of 10 9 ; while the number of IMD events range from 10 4 to 10 6 , depending on the 
energy and polarization of the decaying muon. Then, deficits even smaller than few per 
mille due to the injection of RH neutrinos could be eventually testable for such a large 
statistics. However, per mille deficits may be too small because the highest achievable 
precision in the determination of the flux is expected to be at most of ~ 0.1%. But 
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being conservative, it could be up to a factor ten times larger, and then of the same 
order as the largest possible deficit for the E2 extension. 

At this point one may wonder about the consistency of these large deficits with the 
interpretation of present neutrino oscillation experiments summarized in Table [2j As 



v source 


Experiment 


Detection 


Reactor 
(/? decay) 


ralo Verae |io] . L-riUUZi jlYj 
KamLAND [E] 


~h 


Solar 


SNO [19], Borexino [20] 


Ve -> Vp 


Atmospheric 
(tt and /x decays) 


Super Kamiokande [2"T] 


Vfj, -> I/ T 


Accelerator 
(-7T and K decays) 


K2K [22], MINOS [23], 
CHORUS [21], NOMAD [25] 
MiniBooNE [26] 


^ ~> 


Accelerator 
(fi decay at rest) 


LSND [2] 
KARMEN [27] 


Z^u -> ^ 



Table 2: Neutrino source for the different oscillation experiments and search process. 

we will argue, these seem to be largely insensitive to new four-fermion interactions in- 
volving an electron, a muon and the corresponding neutrinos. Reactor experiments are 
initiated by electron antineutrinos from f3 decay, they are then LH and fully described 
by the SM. Similarly, solar neutrinos have electronic flavour, are also LH and produced 
by SM reactions. Atmospheric neutrinos are decay products of pions and muons from 
cosmic rays, and may include RH neutrinos. However, since the flux of cosmic rays 
is isotropic, and atmospheric neutrino oscillation experiments only compare fluxes of 
muon neutrinos coming from different directions, they are not sensitive to a possible 
deficit in the total number of initial LH neutrinos from muon decays. On the other 
hand, in accelerator experiments looking for — > v T or — > z/ e , neutrino beams are 
mainly formed by muon neutrinos originating from pion decays (with a fairly small 
contamination) [22], and then LH and with the SM oscillation pattern. Finally, in 
accelerator experiments looking for — > V e , the muon antineutrinos are produced in 
/i + decays, and therefore they are sensitive to the NI we are interested in. However, 
what they measure is the number of positrons produced by inverse f3 decay, looking 
for an excess of electron anti-neutrinos instead of looking for a deficit in the observed 
number of muon anti-neutrinos. (The excess reported by the LSND experiment [2] has 
no explanation in this setup, which predicts a deficit.) Hence, there appears not to be 
any contradiction between the significant deficit predicted by the NI considered here 
and the interpretation of current oscillation experiments. 
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5 Further phenomenological implications 



In specific models that contain new fields and interactions there are in general further 
observable effects, as for instance the production of the new particles at large colliders. 
This is the case of the simple model E2 discussed in Section |3j If this type of NI 
saturates the EWPD limits, 

A, 



Vie 



M A [TeV] - ' 4 - (13) 

implying a relatively large X^ e and a light A. The LHC will be able to uncover such a 
scalar triplet for A masses up to 900 GeV and an integrated luminosity of 30 fb _1 [25] 

B 

For this model the relation in Eq. ([7j) translates into a correlation between the scalar 
masses M^a and/or the Yukawa couplings / ee , MM , A^e, namely 



feefw 



M 2 

T) 



32V2G F . (14) 

ivi A 



This implies that for these NI to have sizable effects at a neutrino factory we must have 
relatively large fee,^ and light 77. The main signals in this case are missing energy plus 
one or two leptons £ = e, fi because a charged (neutral) rj has only sizable decays 
into Iv iyv) □. Thus, W, W Z and WW production can provide too large irreducible 
backgrounds for observing these scalars. Other SM processes like ti production can 
also give large backgrounds. The search for this scalar doublet is similar to left slepton 
searches assuming that they only decay into a LH charged lepton and the lightest 
supersymmetric particle [20] • (Although in general sleptons can also have cascade 
decays and be decay products of other supersymmetric particles [3T].) In this case the 
LHC discovery limit is ~ 300 GeV for an integrated luminosity of 30 fb _1 [32], but this 
is assuming that the slepton doublets coupling to the first two families are degenerate. 
However, in our model there is only one scalar doublet, scaling the corresponding limit 
to 250 GeV after correcting by the factor of 2. It may also happen that though the A 
and rj may be too heavy to be directly observed at the LHC (e.g. Ma > 1 TeV and 
> 250 GeV), their effects can still be observable at a neutrino factory provided 
/ and A fulfill Eqs. ([TBI and ffT4"|) : / ee , wt — 2X fie > 0.8. Obviously, large enough 
lepton colliders are better suited for searching these scalars, since they couple mainly 
to leptons, and these colliders allow for a better kinematical reconstruction. On the 
other hand, the relation in Eq. (|14p . as the cancellation of other possible LF violating 
couplings, does not appear to be natural in this simple model, but they could be in 
more complicated frameworks. 



8 This estimate is larger than those quoted in [35J for in this model the scalar triplet only couples to 
e and fi but not to taus, which are more difficult to identify and have larger backgrounds. In contrast 
with type II see-saw models, this triplet is not the only source of neutrino masses, and light neutrino 
masses and mixings can be reproduced even when the A — r conplings vanish. 

9 See [Ml for alternative scalar doublet models. 



11 



Finally, we note that one might consider other SM additions/extensions generating 
Sg^ L , as for instance, heavy neutrino singlets or triplets mixing with the electron or 
muon. In theses cases, however, the relatively small mixings allowed by EWPD [33] do 
not allow contributions large enough to produce a sizable deficit. 
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Appendix 



Our fits are performed using a x 2 analysis for the experimental data collected in Tables [3] 
and HI This data set is essentially the same used in [33], but updating the top [34] 



Quantity 


Experimental Value 


Standard Model 


II It \ V J 


1 7Q 1 _|_ 1 Q 


1 7^ 9 
1 / o.z 


M w [GeV] 


80.399 ± 0.023 


80.373 


2 

91 


0.3005 ± 0.0012 


0.3039 


al 

ilR 


0311 ± 0010 


0.0301 


L 


2.51 ±0.033 


2.46 


Ob. 


4.59 ±0.41 


5.18 


97 


-0.040 ±0.015 


-0.0399 


97 


-0.507 ±0.014 


-0.507 


Qw (Cs) 


-72.74 ± 0.46 


-73.13 


Qw (Tl) 


-116.4 ±3.6 


-116.7 


Qw (e) (M0ller) 


-0.0403 ±0.0053 


-0.0473 


Z^j 1 v ui\ 


1.0000 ±0.0006 


1 


° /°SM 


0.981 ±0.057 


1 



Table 3: Non Z-pole (pseudo-) observables included in the global fits. 



and W [35] mass values as well as the world average for a s [36]. We also include 
light flavour data at the Z-pole [37] , the Tevatron determination of the effective weak 
mixing angle [38], and further low energy observables, such as the weak charges for 
Thalium (atomic parity violation) and electron (M0ller scattering) [TJ, the unitarity 
of the first row of the CKM matrix [39], and obviously, the data from inverse muon 
decay [40]. The SM best fit values for each observable are also included in the Tables. 
They have been computed using ZFITTER 6.43 [4T], complemented with routines for 
the computation of some of the low energy observables. For the minimizations of the 
X 2 we use the package MINUIT [42J. The SM minimum is obtained for the following 
values: 

m H = 96±H GeV, m t = 173.2 ± 1.3 GeV, 

M z = 91.1875 ± 0.0021 GeV, a s (m|) = 0.1184 ± 0.0007, 
AaS (mi) = (276.9 ± 3.3) • 10" 4 . 

The new physics effects are incorporated adding to the SM predictions the tree level 
contributions from the new operators. The impact of this extension on the best fit 
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Quantity Experimental Value Standard Model 



Aa& (Ml) 
as (Mf) 


0.02758 ± 0.00035 
0.1184 ±0.0007 


0.02769 
0.1184 


M z [GeV] 
T z [GeV] 
a° H [nb] 

it L J 

^FB 
^FB 
^FB 


91.1876 ± 0.0021 
2.4952 ± 0.0023 
41.541 ±0.037 
20.804 ±0.050 
20.785 ±0.033 
20.764 ±0.045 
0.0145 ±0.0025 
0.0169 ±0.0013 
0.0188 ±0.0017 


91.1875 

2.4957 

41.479 

20.741 

20.741 

20.788 

0.0164 


A (SLD) 
^ (SLD) 
A T (SLD) 


0.1516 ±0.0021 
0.142 ±0.015 
0.136 ±0.015 


0.1477 


^e(Pr) 
4 T (P T ) 


0.1498 ±0.0049 
0.1439 ± 0.0043 




R° b 

A o,b 

^FB 
aO,c 
^FB 


0.21629 ±0.00066 
0.1721 ±0.0030 
0.0992 ±0.0016 
0.0707 ±0.0035 
0.923 ±0.020 
0.670 ±0.027 


0.2158 
0.1722 
0.1035 
0.074 
0.935 
0.668 


^FB 

A s 

Rs 1 R u +d+s 


0.098 ±0.011 
0.895 ±0.091 
0.371 ±0.022 


0.1037 
0.936 
0.359 


/Ohad 
VFB 

sin 2 ^ 


0.0403 ± 0.0026 
0.2315 ±0.0018 


0.0424 
0.23143 



Table 4: Z-pole (pseudo) observables included in the global fits. 
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values for the SM inputs is insignificant. For the El completion, which allows for the 
largest departure from the SM, 

92^° GeV, m t \ E1 = 173.2 ± 1.3 GeV, 

91.1875 ±0.0021 GeV, a s {m%)\ e1 = 0.1184 ± 0.0007, 

(276.8 ±3.3) ■ 10~ 4 . 

Obviously, this is because none of the new operators directly contributes to (almost) 
any of the considered observables; and they indirectly do only through G M , which is 
in excellent agreement with the SM prediction. Thus, 5gi L tends to be negative in 
order to compensate the effect of the other operators and the SM parameters do not 
feel the presence of the new parameters. Once we have introduced the new operators, 
we perform a scan over the parameter space and reconstruct the probability density 
function (p.d.f.) using an acceptance/rejection method. From the resulting distribution 
we compute the confidence regions at 90% C.L., as well as the marginal p.d.f. for each 
parameter, which we use to obtain the limits and confidence intervals quoted in the 
text. 



El 



M 



z\ei 



El 



15 



References 



C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1; see also 
M. C. Gonzalez-Garcia and M. Maltoni, Phys. Rept. 460 (2008) 1 [arX iv:0704. 18001 
[hep-ph]]. 

A. Aguilar et al. [LS ND Collaboration], Phys. Rev. D 64 (2001) 112007 
|arXiv:hep-ex/0104049l . 

T. Abe et al. [ISS Detector Working Group], JINST 4 (2009) T05001 
[arXiv:0712.4 129 [physics. ins-det]]; for a recent review on muon factories see 

D. M. Kaplan. larXiv:0910.3154l [phvsics.acc-ph], 

W. Buchmuller and D. Wyler, Nucl. Phys. B 268 (1986) 621. 
S. Weinberg, Phys. Rev. Lett. 43 (1979) 1566. 

A. de Gouvea and J. Jenkins, Phys. Rev. D 77 (2008) 013008 [ arXiv:0708. 13441 
[hep-ph]]. 

F. del Aguila, S. Bar-Shalom, A. Soni and J. Wudka, Phys. Lett. B 670 (2009) 
399 [arXiv:0806.0876l [hep-ph]]. 

A. Aparici, K. Kim, A. Santamaria and J. Wudka, Phys. Rev. D 80 (2009) 013010 
|arXiv:0904.3244l [hep-ph]]. 

P. Langacker and D. London, Phys. Rev. D 39 (1989) 266. 

M. B. Gavela, D. Hernandez, T. Ota and W. Winter, Phys. Rev. D 79 (2009) 
013007 |arXiv:0809.345T1 [hep-ph]]. 

W. Fetscher, H. J. Gerber and K. F. Johnson, Comparison with Phys. Lett. B 173 
(1986) 102; C. A. Gagliardi, R. E. Tri bble and N. J. Williams, Phys. Rev. D 72, 
073002 (2005) [arXiv:hep-ph/0509069| . 

C. Biggio, M. Blennow and E. Fernandez-Martinez, JHEP 0908 (2009) 090 
[arXiv:0907.0097l [hep-ph]]. 

G. Prezeau and A . Kurylov, Phys. Rev. Lett. 95 (2005) 101802 
[arXiv:hep-ph/0409193] . 



E. Ma, Phys. Rev. Lett. 86 (2001) 2502 |arXiv:hep-ph/0011121| . 



A. Abada, C. Biggio, F. Bonnet, M. B. Gavela and T. Hambye, JHEP 0712 
(2007) 061 [ arXiv:0707.4~058l [hep-ph]]; F. del Aguila, J. A. Aguilar-Saavedra, 
J. de Bias and M. Zralek, Acta Phys. Polon. B 38 (2007) 3339 |arXiv:0710.2"923l 
[hep-ph]]; F. del Aguila, J. A. Aguilar-Saavedra, J. de Bias and M. Perez- 
Victoria, arXiv:0806.1023 [hep-ph]; see also T. Ohlsson, T. Schwetz and H. Zhang, 
larXiv:0909.0455l [hep-ph]. 



16 



F. Boehm et al, Phys. Rev. D 64 (2001) 112001 |arXiv:hep-ex/0107 009 1 . 

M. Apollonio et al. [C HOOZ Collaboration], Eur. Phys. J. C 27 (2003) 331 
[arXiv:hep-ex/0301017] . 

S. Abe et al. [KamLAND Collaboration], Phys. Rev. Lett. 100 (2008) 221803 
[arXiv: 0801 .45891 [hep-ex] ] . 

B. Aharmim et al. [S NO Collaboration], Phys. Rev. C 72 (2005) 055502 
[arXiv:nucl-ex/0502'02ll . 

C. Arpesella et al. [Borexino Collaboration]. Phys. Lett, B 658 (2008) 101 
[arXiv: 0708 .2251 [astro-ph]]. 

Y. Ashie et al. [Super-K amiokande Collaboration], Phys. Rev. D 71 (2005) 112005 
[arXiv:hep-ex/0501064] . 

M. H. Ahn et al. [K 2K Collaboration], Phys. Rev. D 74 (2006) 072003 
|arXiv:hep-ex/0606032] . 

D. G. Michael et al. [M INOS Collaboration], Phys. Rev. Lett. 97 (2006) 191801 
[arXiv:hep-ex/0607088] . 

E. Eskut et al. [CHORUS Collaboration], Nucl. Phys. B 793 (2008) 326 
[arXiv:0710.336i1 [hep-ex] ] . 

P. Astier et al. [NO MAD Collaboration], Nucl. Phys. B 611 (2001) 3 
[arXiv:hep-ex/0106102] . 

A. A. Aguilar-Arevalo et al. [The MiniBooNE Collaboration], Phys. Rev. Lett. 98 
(2007) 231801 [arXiv:0704.1500l [hep-ex]]. 

B. Armbruster et al. [K ARMEN Collaboration], Phys. Rev. D 65 (2002) 112001 
[arXiv:hep-ex/020302l] . 

F. del Aguila and J. A. Aguilar-Saavedra, Nucl. Phys. B 813 (2009) 22 
[arXiv: 08 08 .2468 [hep-ph]]; see also K. Huitu, J. Maalampi, A. Pietila and 
M. Raidal, Nucl. Phys. B487, 27 (1997) [ |hep-ph/960631ip ; J. F. Gunion, 

C. Loomis and K. T. Pitts, [ ]hep-ph/96 102370 ; A. G. Akeroyd and M. Aoki, Phys. 
Rev. D72, 035011 (2005) Qhep-ph/0506176P ; A. Hektor, M. Kadastik, M. Mun- 
tel, M. Raidal and L. Rebane, Nucl. Phys. B787, 198 (2007) HarXiv: 0705. 14951 
[hep-ph]]; P. Fileviez Perez, T. Han, G. Y. Huang, T. Li and K. Wang, 
Phys. Rev. D78, 071301 (2008) HarXiv: 0803 . 34501 [hep-ph]]; P. Fileviez Perez, 
T. Han, G. y. Huang, T. Li and K. Wang, Phys. Rev. D78 015018 (2008) 
HarXiv: 0805 .35361 [hep-ph]]. 

[29] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 102 (2009) 041801 
[arXiv: 0809. 4903 [hep-ex]]; and references there in. 



17 



[30] F. del Aguila and L. Ametller, Phys. Lett. B 261 (1991) 326. 

[31] H. Baer, C. h. Chen F. Paige and X. Tata, Phys. Rev. D 49 (1994) 3283 
|arXiv:hep-ph/9311248]; H. Baer, B . W. Harris and M. H. Reno, Phys. Rev. D 57 
(1998) 5871 |arXiv:hep-ph/9712315] . 

[32] Yu. M. Andreev, S. I. Bityuk ov and N. V. Krasnikov, Phys. Atom. Nucl. 68 (2005) 
340 [Yad. Fiz. 68 (2005) 366] [arXiv:hep-ph/0402229| ; see for a review F. del Aguila 
et al, Eur. Phys. J. C 57 (2008) 183 |arXiv:0801. 18001 [hep-ph]]. 

[33] F. del Aguila, J. de Bias and M. Perez- Victoria, Phys. Rev. D 78 (2008) 013010 
jarXiv: 0803 .40081 [hep-ph]]; F. del Aguila, J. A. Aguilar-Saavedra and J. de Bias, 
larXiv:0910.2720l [hep-ph]. 

[34] [Tevatron Electroweak Working Group and CDF Collaboration and DO Collab], 
larXiv:0903.2503l [hep-ex]. 

[35] [Tevatron Electroweak Working Group and CDF Collaboration and DO Collab], 
larXiv:0908. 13741 [hep-ex]. 

[36] S. Bethke. la?Xiv:0908. 11351 [hep-ph], 

[37] [ALEPH Collaboration and DELPHI Collaboration and L3 Collaboration and ], 
Phys. Rept. 427 (2006) 257 |arXiv:hep-ex/0509008] . 

[38] D. E. Acosta et al [CDF Collaboration], Phys. Rev. D 71 (2005) 052002 
[arXiv:hep-ex/041105 9|; V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 
101 (2008) 191801 |arXiv:0804.3220l [hep-ex]]. 

[39] J. C. Hardy and I. S. Towner, Phys. Rev. C 79 (2009) 055502 [ arXiv:0812. 12021 
[nucl-ex]]. 

[40] S. R. Mishra et al, Phys. Lett. B 252 (1990) 170. 

[41] A. B. Arbuzov et al, "ZFITTER: A semi-analytical program for fermion 
pair production in e+ e- Comput. Phys. Commun. 174, 728 (2006) 
|arXiv:hep-ph/0507146] . 

[42] F. James and M. Roos, Comput. Phys. Commun. 10 (1975) 343. 



18 



